A Solvable Model for Homopolymers and the Critical 

Phenomena 



M. Cranston* L. Koralovj S. Molchanovf B. Vainberg^ 



Abstract 

We consider a model for the distribution of a long homopolymer with a zero- 
range potential at the origin in R 3 . The distribution can be obtained as a limit of 
Gibbs distributions corresponding to properly normalized potentials concentrated 
in small neighborhoods of the origin as the size of the neighborhoods tends to zero. 
The distribution depends on the length T of the polymer and a parameter 7 that 
corresponds, roughly speaking, to the difference between the inverse temperature in 
our model and the critical value of the inverse temperature. 

At the critical point 7 cr = the transition occurs from the globular phase 
(positive recurrent behavior of the polymer, 7 > 0) to the extended phase (Brownian 
type behavior, 7 < 0). The main result of the paper is a detailed analysis of the 
behavior of the polymer when 7 is near 7 cr . 

Our approach is based on analyzing the semigroups generated by the self-adjoint 
extensions £ 7 of the Laplacian on Cq°(M 3 \{0}) parametrized by 7, which are related 
to the distribution of the polymer. The main technical tool of the paper is the 
explicit formula for the resolvent of the operator £ 7 . 

Key words: Gibbs measure, homopolymer, zero-range potential, phase transition, glob- 
ular phase, diffusive phase. 
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1 Introduction 

The study of polymer models has been a very active area of research in mathematical 
physics in recent years. Many of the physically relevant problems have been outlined in the 
paper of Lifschitz, Grosberg and Khokhlov [17) . To name just a few recent articles on the 
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subject, with apologies to the many authors who have been omitted, we cite [2], [5], [I], [6], 
[7J, [8]. In addition, there is an interesting exposition with many valuable references in [12]. 
The approach to the subject involves many essential ideas from statistical mechanics 
including strong connections to the developments in [UJ, [13] and [IB] . 

Let us give a general qualitative description of the problem. We start with a path 
space Q = C([0, T], IR 3 ), a Hamiltonian if^ : f2 — > (— oo,0] and a probability measure 
on fi. A continuous function w G O will be thought of as a realization of the polymer. 
The parameter t G [0,T] can be intuitively understood as the length along the polymer 
(although the functions ui = u(t) are not differentiable and the genuine notion of length 
can not be defined). One can consider models with the Hamiltonian given by 



cT 



H T (u) — — / / v(u(si) — u(s 2 ))ds 1 ds 2 , 
Jo Jo 



where v : IR 3 x IR 3 — ► [0, oo) is a local attractive potential (see [3], [H]). Our paper 
concerns a simpler "mean field" type model, where the polymer chain interacts with the 
external attractive potential (as in p2]). Namely, consider the Hamiltonian 

H T (cu) = - [ v(cu(s))ds, (1) 



where v is a nonnegative compactly supported potential that is not identically equal to 
zero. The model considered here is that of a homopolymer since the potential is time- 
independent. Other models (of heteropolymers) which we do not consider in this paper 
allow for a potential that is a random stationary function of time (see [2] for example). 

Let Pf. be the Wiener measure on Q = C([0,T],IR 3 ) shifted so that the trajectories 
start at x almost surely. This will be the reference measure corresponding to the infinite 
temperature, i.e., to the inverse temperature j3 = 0. For a value of inverse temperature 
(3 > 0, the polymer is distributed according to the Gibbs measure P^ T , whose density 
with respect to Pf. is 

The normalizing factor Zp(T,x) = E^e~ /3HT is called the partition function. 

The phase transitions for the polymers are well-understood at the physical level (see 
[TT]). For large (3 (/3 > (3 cr for some (3 cr > 0) and any t G [0, T], with high probability 
with respect to the measure P% T the values of u(t) are bounded by a constant that is 
independent of t. Thus the polymer is in the globular state. However, for f3 < f3 cr the 
typical shape of oo(t), t G [0, T], is that of a Brownian path in M 3 , and the polymer is 
said to be in the diffusive state. Rigorous results of this nature have also been proven 
in some generality (see [12], [2] for example). Of particular interest is the behavior of 
the polymer for the values of (3 that are near (3 cr . Unfortunately, the detailed analysis of 
critical phenomena for general Hamiltonians is largely outside the range of techniques of 
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modern mathematical physics. Most of the results on this type of problems concern the 
existence of a phase transition. 

However, in our situation (homopolymer with a compactly supported non-negative 
potential), the complete analysis of the critical phenomena is possible. In [9] we studied 
the prevalent behavior of the polymer with respect to the measure P^ T as T — ► +00. In 
particular, we saw that for d > 3 there is f3 cr > such that: 

(a) For (3 > j3 cr , each constant s > and each function S(T) such that S(T) — ► +00 
and T — S(T) — > +00 as T — ► +00, the distribution of uj(S(T) +t), t £ [0, s], with respect 
to Pp T , converges to a recurrent Markov process on [0, s]. 

(b) For (3 < (3 cr , the distribution of u>(tT)/yT, t £ [0, 1], converges to the Brownian 
motion on the interval [0, 1]. 

(c) For (3 = (3 cr) the distribution of u>(tT)/y/T, t £ [0, 1], converges to the distribution 
of a non-Gaussian Markov process on the interval [0, 1]. 

Some of the earlier results for the discrete model based on the random walk and 
zero-range potential can be found in [TO] . 

In the current paper we study the distributions of the continuous polymer in M 3 
which correspond to Schrodinger operators with zero-range potential. The potential is, 
in some sense, supported at the origin. Note, however, that the Hamiltonian H?, the 
corresponding Gibbs measure and the Schrodinger operator Lp defined by ([I]), (T5]) and 
([3]) are meaningless if v is a delta-function concentrated at the origin. Before we define 
these objects, let us note that in the case of a smooth potential, the finite-dimensional 
distributions of T can be expressed in terms of the fundamental solutions ppit, x, y) of 
the parabolic equation 

— = Lpu := -Au + (3vu. (3) 
Indeed, by (T5J) and the Feynman-Kac formula, 

Zp(T,x) = / p p (T,x,y)dy, 
Jr 3 

and 

P^ T (w(ti) £ A u w(t k ) £ A k ) = (4) 
Zn l {T,x) I ... j \ pp(ti,x,xx)...pp(t k - t k _i,x k -i,x k )pp{T - t k ,x k ,y)dydx k ...dx 1 , 

J Ai JA k Jr3 

where k > 1, < t\ < ... < t k < T and Ai, A k are Borel sets in R 3 . Here we use the 
convention that j9g(0, x, y) = 5 x (y). The measure P^ T corresponds to a non- homogeneous 
Markov process on [0,T] (see [9]). 

Now, instead of the operators Lp, we start with the Schrodinger operators with zero- 
range potential. They are defined (see [TJ) as self-adjoint extension in L 2 (IR 3 ) of the 
operator 

A/2 : C™(R 3 \ {0}) -f C™(R 3 \ {0}). 
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There is a family of such extensions £ 7 which depend on a parameter 7 6 R that plays 
the same role as the difference between j3 and (5 cr in the case of smooth potentials (see 
Section[2]). Then we replace operator Lp in ([3]) by the operator £ 7 . This allows us to define 
the finite dimensional distributions and construct the Gibbs measure which corresponds 
to the zero-range potentials. Namely, let p 7 (t, x, y), t > 0, be the kernels of the operators 
in the semi-group generated by £ 7 . We can now define 

Z>y(T,x) = / p {T,x,y)dy, 

and use an analogue of (j3J with p 7 instead of pp and Z^(T,x) instead of Zp(T,x) in 
order to define the finite-dimensional distributions of the measure P 7 t- (While neither 
p (ti,x, Xi) nor Zj(T,x) are defined when x = 0, we can make sense of the expression 

X (T, 0)p 7 (ti, 0, %i) by taking the limit of Z^ X (T, x)pJti, x, X\) as |x| { 0.) It is not 
difficult to check that there is a measure on C([0, T], M 3 ) with such finite-dimensional 
distributions, and that this measure defines a non- homogeneous Markov process on [0, T]. 

There is another way to obtain the measures P 7 T corresponding to zero-range poten- 
tials. One can construct P 7)T as the weak limit of measures corresponding to bounded 
compactly supported potentials as their supports shrink and the values of the potentials 
increase in a particular way. Namely, it turns out that £ 7 , 7 e 1, can be obtained as a 
strong resolvent limit, as e [ 0, of Schrodinger operators £ 7 = A/2+t> 7 with the potentials 
(seeffl) 

f 1 71-2 7\ i X \ 11 11 47T 

«7 = + 7M7). IHIW) = y. (5) 

It will be shown in this paper that the resolvent convergence of the Schrodinger operators 
implies the convergence of the corresponding Gibbs measures. Consider the Hamiltonian 
H^ t given by (p]) with t> 7 instead of v, and the Gibbs measure P 7 ^ given by (|2J) with 
(3 = 1 and T instead of Ht- For each 7 e R and T > there are limits 

understood in the sense of weak convergence of measures on C([0, T], M 3 ) (see Section [2]). 

Unlike the case of a generic smooth potential, all the relevant analytic quantities can 
be found explicitly in the case of the zero-range potential polymer model. This allows us 
not only to obtain an analogue of the results (a)-(c) above, but also provide a detailed 
analysis of the behavior of the polymer distribution when 7 tends to its critical value and, 
simultaneously, T — > +00. 

We recall the main facts about the self-adjoint extensions and prove the convergence 
of P 7 '^ to P T in Section [2j We refer the reader to [lj for a detailed treatment of the self- 
adjoint extensions of the Laplacian on C^°(R 3 \ {0}). Here we only mention that for d > 4 
the only closed self-adjoint extensions of the Laplacian on C^°(M. d \ {0}) is the Laplacian 
on the entire space, and thus there is no analog of the distribution corresponding to a 
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zero-range potential in d > 4. For d = 1,2, the corresponding homopolymer models do 
not exhibit a phase transition, so we do not treat them here. 

The goal of the current paper is to analyze the behavior of the polymer under the 
measure P T for various values of 7 and large T. Of particular interest is the behavior 
of the polymer when 7 is near zero. We provide a detailed analysis of the behavior of 
P 7 (tyt when 7 = 7(T) — ► as T — > +00. 

Two qualitatively different cases can be considered. The analysis of both cases relies 
on a self-similarity property of the measures P T . In the first case, j(T) is bounded 
and such that 7(T)v / T — > +00 as T — ► +00. We shall see that for each constant s > 
and each function S(T) such that j(T) S(T) — > +00 and 7(T)a/T — S(T) — > +00 as 
T — > +00, the distribution of 7(T)a;(S , (T) + (j(T))~H), t G [0, s], with respect to P 7 (t),t> 
converges to a recurrent Markov process on [0, s]. The generator of the Markov process 
will be written out explicitly (see Section [3]). In particular, the radial part is a diffusion 
process on the positive semi-axis with reflection at the origin. These results are also 
applicable in the situation when 7 > does not depend on T. 

In the second case 7(T) is such that 7(T)v / T — > x G [—00, +00) as T — > +00. We 
will show that the distribution of u(tT)/\/T : t G [0, 1], converges to the limiting measure 
on C([0, 1], M 3 ). The limiting measure can be identified as P^i if x > — 00 and the 
Wiener measure on C([0, 1], IR 3 ) if x = —00. In particular, the distribution of uj(tT)/y/T, 
t G [0, 1], converges to the Wiener measure on C([0, 1], M 3 ) if 7 < does not depend on T. 
The limiting distribution for u(T)/\/T can be written out explicitly - it turns out to be 
compound Gaussian (see Section Hj). 

The paper is organized as follows. In Section [2] we recall some facts about the self- 
adjoint extensions of the Laplacian and introduce the associated family of processes with 
the zero-range potential. In Sections [3] and H] we prove the main results concerning the 
cases when j(T)\/T — > +00 and 7(T)a/T — > x G [—00, +00) as T — > +00, respectively. 
Finally, in Section Owe prove the convergence of the processes with potentials t> 7 to those 
with a zero-range potential and prove some technical lemmas regarding the tightness of 
certain families of processes. 

2 Measures Corresponding to Zero- Range Potentials 

All the self-adjoint extensions of A/2 are described in the following theorem, whose proof 
can be found in [I]. 

Theorem 2.1. All the self-adjoint extensions of the Laplacian acting on C^(M. 3 \ {0}) 
to an operator acting on L 2 (M. 3 ) form a one-parameter family £ 7 , 7 G R. The spectrum 
of L 1 is given by 




7 > 



spec(£ 7 ) = (—oo,0], 7 < 0. 
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The kernel of the resolvent of £ 7 is given by 

e -V2\\x-y\ i e -V2X(\x\ + \y\) 

R\,j(x, y) = — : + -= o i 1 1 i » A i s P ec (£ 7 )- 

7r|s-y| V2A-7 27r|x||j/| 

7/7 > 0, i/ien 7 2 /2 zs a simple eigenvalue of with the eigenfunction 

^7 e" 7 ^! 



^ 7 (ar) 



7 27r \x\ 



Since the spectrum of £ 7 is bounded from above, the operators exp(t£ 7 ), t > 0, are 

7) 



bounded in L 2 (R 3 ). The kernel of exp(t£ 7 ), t > 0, is given by 



If e -\x-y\ 2 /2t l r e -V2\(\x\+\y\)+\t 

p(t,x,y) = — e t R Xn (x ) y)d\ = — — — + — - / — = ——d\, (6) 

2tt? y r(a) (27rt)^ 2 4tH2 y r(a) (v^A - 7) |x| \y\ 

where x,y 7^ 0, a > 7 2 /2 and T(a) is the contour in the complex plane that is parallel to 
the imaginary axis and passes through a. Thus pJt,x,y) can be interpreted as a formal 
fundamental solution of the equation 

du 

- = C,u. 

Let 

Z 7 (t,x) = / p (t,x,y)dy, 
where t > 0, x 7^ 0. Formula ([6]) implies that 

Z 7 (t, x) = 1 + — / e A '^= ^^A. (7) 

7 2Wr(«) v^-7 A M 

Define the measures P^ T , x G M 3 , via their finite-dimensional distributions 

P^ jT (w(ti) G Ai, ...,uj(t k ) G A k ) = (8) 



where > 1, < ii < ... < tk < T and A\, ...,A^ are Borel sets in M 3 . Here we use the con- 
ventions that p 7 (0, x, y) = 8 x (y) and Z y (T, 0)p 7 (£i, 0, x x ) = lim^^o ^ 7 (T, x)p 7 (*i, x, xi). 
It is not difficult to show that there indeed is a measure on C([0,T],1R 3 ) with the finite- 
dimensional distributions given by (jHJ). We don't prove this here since the same conclusion 
follows from the proof of Theorem 12.31 below. Note the following self-similarity property. 
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Theorem 2.2. For a > 0, let f a : C([0,T],R 3 ) -> C([0, aT], R 3 ) map a function uj to 
the function f a uo via (f a u)(t) = y/au(t/a). Let / q P 7 t ^ e ^ e push-forward of the measure 
rp by this mapping. Then 

p 7 (to, xy/a, yy/a) = a~ 3/2 p jVE (t, x, y), Z^ta, xy/a) = Z^(t, x) (9) 

and 

faK T = p;f^ Ta . (io) 

Proof. The first statement follows after making the change of variables A' = A/a in the 
integrals in the right hand sides of (EJ) and (J7J) . The second statement now follows from (JHJ) . 

□ 

Let us discuss an alternative construction of the measure PZ't that uses short-range 
potentials. Let v e be defined by (JSJ). Let the Hamiltonian H^ T be given by ([T]) with 
instead of v, and the Gibbs measure be given by (j2J) with (3 — 1 and H^ T instead 
of Ht- The following theorem will be proved in Section 

Theorem 2.3. Let x G R 3 \ 0, 7 6 I and T > 0. The measures P*^ converge, in the 
sense of weak convergence of measures on C([0, T], M 3 ) ; as £ J, 0, to £/ie measure P^t- 

j.' 

From (|S]) it follows that the finite-dimensional distributions of P / T converge to those 

of P° T as x' — > and 7' — > 7 G JR.. If 7' — ► —00, then they converge to the finite- 
dimensional distributions of the Wiener measure on C([0, T], K 3 ). From Theorem 12.31 and 

x ' 

Lemma EH] (from Section E]) it follows that the family {Py T , 7' < c, < 1} is tight for 
each c6R, and therefore 

lim P~:; T = PV lim P X ' T = P° T , (11) 

where P^ is the Wiener measure and the limits are understood in the sense of weak 
convergence of measures. 



3 Distribution Above the Critical Point 

In this section we assume that 7 = j(T) is bounded and such that j(T)\/T — > +00 
as T — > +00. In particular, this covers the case when 7 > does not depend on T. We 
examine the behavior of the polymer paths with respect to P 7 m T as T — > +00. First, 
we need the asymptotics of Pi(t, x, y) and Zi(t,x) when t — > +00. 

Lemma 3.1. We have the following asymptotic expressions: 

Pi(t,x,y) = ex.p(t/2)(ip 1 (x)ip 1 (y) +q(t,x,y)), (12) 
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where lim 4 ^ +00 supui i y i >fe \q(t,x,y)\ = for each k > 0; 

Zxit.x) = exp(t/2)||^i|| L i( M 3)(^i(rr) + Q(t,x)), (13) 
where lim 4 ^ +00 supu^j. \Q(t, x)\ — /or each k > 0; 

{s,x)p 1 {t,x,y) = exp((t- s)/2)\\ij; 1 \\l^ R3) (ij; 1 (y) + q(s,t,x,y)), (14) 

w/iere lim Sji ^ +00 sup|as|<fc-i,| v |>fc l?( s > *> x > v)\ = / or eac/i k > °- 

Proof. Formula (fT2"j) follows immediately from (jH]) if the integral over T(a) in is re- 
placed by the integral over T(b), < b < 7 2 /2, plus the residue at A = 7 2 /2 = 1/2. 
Formula (fT2"|) follows similarly from (JTj). □ 

Next, let us study the distribution of the end of the polymer with respect to the 
measure P 7 ( T ) jT as T — > +oo. 

Theorem 3.2. // 7(T) zs bounded and such that 7(T)v / T — > +oo as T — > +oo ; i/ien 
i/ie distribution of j(T)uj(T) with respect to the measure P 7 ( T ) )T converges, weakly, as 
T — > +oo ; to £/ie distribution with the density V'l/IIV'ilU 1 ® 3 )- 

Proof. For fixed the density of a;(T) with respect to the Lebesgue measure is equal to 

z 7 (t)( t , x)p l(T) (T,x,y). 
Therefore, the density of j(T)u(T) is equal to 

7(T)~ 3 Z^ ( 1 T )(T, x)p 7 ( T )(T, x, y/^{T)) = Z 1 ^ 1 (7 2 (T)T, x7(T))p 1 (7 2 (T)T, ary(T), y). 
Here we used (jUJ) with 7 = 1, a = 7 2 (T) and t = T. The latter expression is equal to 

I |Vi IIZiV)^ 1 ^) + 9(T 2 (r)r, 7 2 (T)T, S7(T), y)), (15) 



where g is the same as in Lemma [331 When T — > +00, the expression in (fT5l) converges to 
ipi(x) /\\iPi\\li(m,3} uniformly in \y\ > k by Lemma I3TT1 This justifies the weak convergence. 

□ 

Now let us examine the behavior of the polymer in a region separated both from zero 
and T. Let S(T) be such that 

lim 7(T)JS(T) = lim 1 {T)JT - S(T) = +00. (16) 

Let s > be fixed. Consider the process y T {t) = ^(T)u{S{T) + t/7 2 (T)), < t < s. Let 

/. \ Pi(t,x,y)ipi(y) , , 
r(t,x,y) = —— exp(-i/2), y ^ 0. (17) 

Observe that / R3 r(t, x, y)dy = 1, and therefore r serves as the transition density for a 
Markov process. Also note that ipf serves as an invariant density for the process. 
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Theorem 3.3. If ^(T) is bounded and such that r ){T)\/T — > +00 as T — > +00 ; then the 
distribution of the process y T (t) with respect to the measure P 7 m,T converges as T — > +00, 
weakly in the space C([0, s], M 3 ), to the distribution of a stationary Markov process with 
transition density r(t,x,y) and the invariant measure whose density is 

Proof. Let < t 1 < ... < t n < s. The density of the random vector (y T (ti), y T {t n )) 
with respect to the Lebesgue measure on M 3n is equal to 

p T (x 1 ,...,x n ) = 

7(T)- 3 "(I 7(r) (T,x))-X (r) (S(T) + ^)^7(T)(^y> ^ry ^j) - 

_ ,tn t n —i X n _i X n — . c( r r\ ^ n \ 

7v ; 7 2^ T ) j(T) 7(T) 7 2 (T) 7 (T) 

Applying (jUj) with 7 = 1, a = 7 2 (T) and t = T, we see that the right hand side of this 
equality is equal to 

(Z 1 ( 7 2 (T)T, a;7(T)))- 1 p 1 (7 2 (T)5(T) + h, ^{T^p^h - h, x x , x 2 ) ... 

• •• Pxitn - t n _i,a; n _i,x„)Zi(7 2 (T)(T - S(T)) - t n ,x n ). 

We replace here all factors p 1( except the first one, by r using (TTTT) . We replace the factors 
(Z 1 (7 2 (T)T,x7(T)))- 1 p 1 (7 2 (T)5(T)+t 1 ,x 7 (T),x 1 )andZ 1 (7 2 (T)(T-5(T))-t n ,x„) by 
their asymptotic expansions given in Lemma 13.11 This leads to 

p T (x 1 , x n ) = ipl(xi)r(t 2 - h, x 1 ,x 2 )...r(t n - £ n _ 1; x n -i,x n ) + o(l), T -> +00, 

where the remainder tends to zero uniformly in (xi,...,x n ) with min(|xi|, |x n |) > k. 
Since > can be chosen to be arbitrarily small, this justifies the convergence of the finite 
dimensional distributions of y T to those of the Markov process. It remains to note that 
the family of measures induced by the processes y T is tight, as follows from Lemma [531 □ 

Remark. Let us describe the generator of the limiting Markov process. Namely, let 
C*o(M 3 ) be the space of continuous functions on M. 3 with a finite limit at infinity. Consider 
the differential operator 

> r / % 1a , \ (Vt/'ifx), Vu(x)) 1 . , . 1 .du(x) 

Mu{x) = -Au(x) + V ^ U ; r V ;; = -Au{x) - 1 + i-r)-^, x ^ 0, 
2 wi\ x ) 2 \x\ or 

Mu{0) = lim(^Au(x) - (1 + _L) 9u Wn 



\x\10 2 \x\ dr 

Note that in the spherical coordinates M can be written as 

1 ,<9 2 w 1 ~ du 

2 Q r 2 r 2 q t ' 
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where A is the Beltrami-Laplace operator on the sphere. For fixed x, the function r(t, x, y) 
satisfies 

9r{t £ y) =M*r(t,x,y), t>0 )2/ 6R 3 \0, 

where M* is the formal adjoint of M. From here it easily follows ([15]) that the generator 
in the space C*o(M 3 ) of the Markov family with transition density r(t, x, y) is given by the 
operator M with the domain 

V = {u G C 2 (M 3 \ {0}) n C (R 3 ) : Mu G C (M 3 ), lim ^ M(r ' y) ^ (y) ~ m(0) = 0}, 

where /x is the Haar probability measure on S 12 . The radial part of the limiting process is 
then a diffusion on the positive semi- axis with unit diffusion coefficient, unit drift towards 
the origin and reflection at the origin. If we denote the radial part of the limiting process 
by Rt, then the spherical part of the limiting process St on S 2 satisfies 

dS t =—dB t , R t ^0. (18) 
tit 

Here B t is the diffusion on the sphere whose generator is one half of the Beltrami-Laplace 
operator. It is not difficult to see that if we denote the the first time when R t = by r, 
then the set of limit points of Sf as t f r coincides with the entire sphere S 2 . For t > 
the distribution of S T+t is uniform on the sphere. 

Remark. If instead of assuming that liniT->+oo l(T) a/ S(T) = +oo we assume that 
S(T) = 0, the result of Theorem 13.31 will hold with the only difference that the initial 
distribution for the limiting Markov process will now be concentrated at x liniT^+oo j(T), 
instead of being the invariant distribution, provided that the latter limit exists. 



4 Distribution Near and Below the Critical Point 

In this section we assume that 7 = 7(T) is such that ^(T)VT — > k G [—00, +00) as T — > 
+00. In particular, this covers the case when 7 < does not depend on T. We examine 
the behavior of the polymer paths with respect to P 7 m t as T — > +00. 

Consider the process y T (t) = u(tT)/ VT, < t < 1. Applying ffTUl) with a = 1/T, we 
see that the distribution of the process y T with respect to the measure P 7 m,T coincides 

with the distribution of the process u(t), t G [0, 1], with respect to the measure P 7 ( T ) N /y r 

Since x/y/T — > and "f(T)\/T — > x as T — > +00, the following theorem is a consequence 
of dHJ>. 

Theorem 4.1. J/ 7 = 7(7") is swc/i i/iai 7(T)V / T — > x G (—00, +00) as T — > +oo ; 
i/ien i/ie distribution of the process y T (t) with respect to the measure P 7 ( T ) T converges as 

T — > +00 to £/ie measure p" , . 
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If 7 = is such that r )(T)\/T — > — oo as T — > +00, i/ien i/ie distribution of the 

process y T {t) with respect to the measure P 7T converges as T — > +00 to i/ie distribution 
of the 3- dimensional Brownian motion. 

The finite-dimensional distributions of the limiting process can be written out using 
(J6j) and (jBJ). In particular, when x = 0, the distribution of the limiting process is the 
same as the distribution of the corresponding process in the case of the smooth potential 
with (3 = f3 cr (see [9]). 

Let us show that the limiting distribution of y T (T) is compound Gaussian. 

Theorem 4.2. The distribution ofu(l) (induced by the measure x ) is compound Gaus- 
sian, i.e., its density is given by 

qM = l (2^F^ (r)dr ' 

where v„{t) > 0, r E [0, 1], x G R. 

Proof. By (fTTj) . ((Tj) and (|HJ), the density of u;(l) with respect to the Lebesgue measure is 
equal to 



q x (y) = hm -=— — - = — / — = — dX 



Nio Z„(l,x) 2tt J r(a) (y/2X- H)\y\ \Jv{a) (V2A - x)A 
where a > x 2 /2. Note that 

-e~ Xr dr = — — -—dr, X G Via). 



2n\y\ J (2ttt) 3 / 2 " J dr \(2ttt) 3 / 2 I X 

We put this expression in the first integral in the previous formula. The integration by 
parts above was needed to guarantee the absolute convergence of the double integral and 
to change the order of integration. Below we will move the derivative back to the second 
factor. Thus 



00 



d I e -\y?/iT \ r e A(1 ~ r) (iA 

q*(y) = I -r 77, — uTi u„(T)dT, u„(t) = K / — = — , (19) 

Jo dr y(27rT) 6 ^J J r(a) (V2A - x)A 

where 

e x dX 



K 

'r(a) (V2A - x)A, 

If r > 1 and ReA > a, then the integrand in the second integral in (1191) is analytic in A 
and of order 0(|A|~ 3 / 2 ) as |A| — > 00. Thus u h (t) = for r > 1, and 

«M = {j2^ ) U ^ T)dT = I (2^^ (T)rfT ' (20) 
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where 



v«(t) = -u'Jf) = K / — = -, < r < 1. 

Jr(a) (V 2A - x) 

Since the integral over T(a) is equal to the integral over a contour around negative A- 
semiaxis (plus the contribution from the residue at x 2 /2 if x > 0), we have 



2i / + 527iixe^ {1 - T)/2 , < r < 1, (21) 



r(a) (V2X - x) 7o 2a + x 2 

where 5 = 1 if x > 0, 5 = if x < 0. 
Let us evaluate K. We have 



x _ r e x d\ _ r {V2X + x)e x dX f y/2Xe x dX C xe x dX 

k = J m (V2X-h)X ~ Jr(a) (2A - x 2 )A = J r[a) (2A - x 2 )A + j T[a) (2A-x 2 )A' 

(22) 

The first integral Ji in the right-hand side can be written as the integral over the contour 
around the negative A-semiaxis plus the contribution from the residue at x 2 /2, i.e. 



oo 



h = -2V2i / — — + 2m—-. 

Jo (2o- + x 2 ) v /o" \x\ 

We make the substitution a — *> a 2 in the first term of the right-hand side and rewrite the 
second term using the identity 

7T f°° y/2da 



|x| J (j 2 + x 2 /2 

The second integral in the right-hand side of (j22p can be expressed through the residues 
at A = and A = x 2 /2. This implies, that 

K- 1 = 2^/2i\ — v —da + 2m . 

J a 2 + x 2 /2 x 

This and (|2"Tj) show that vj^r) > for all < r < 1, — 00 < x < 00. This together with 
f l20|) justify the statement. □ 



5 Proof of Theorem 12. 3L Tightness of Certain Fami- 
lies of Processes 

Proof of Theorem \2.3[ As follows from ([3]) and (jl]), the finite-dimensional distributions of 
P^y are given by 

P%.(w{t 1 )eA 1 ,...,u(t k )eA k ) = 
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(Z*) 1 (T,x) / ... / / ^(ti,x,xi)...p^(t fc - tk-i,x k -i,Xk)^f(T - t k ,x k ,y)dydx k ...dx 1 , 

J A! J A k JR 3 

where pt(t, x, y) is the fundamental solution of the parabolic equation 

du ~ F 1 . 

and 

Zl(t,x) = / pZ(t,x,y)dy. 

Since P T (oj(ti) = for some 1 < i < k) = 0, in order to demonstrate the conver- 
gence of the finite-dimensional distributions, it is sufficient to prove that lim e ^ Z^(t, x) = 
Zj(t,x) and that for continuous functions with compact supports (f\,...,(p k such that 
supp(y?j) cR 3 \0, 1 < i < k, we have 



lim / ... / ff(t ll x,x 1 )...p £ Jt k -t k _ 1 ,x k ^ 1 ,x k )(pi(x 1 )...(p k (x k )dx k ...dx 1 = (23) 
/ ... / p^(t 1 ,x,x 1 )...p 1 (t k -t k ^ 1 ,x k ^ 1 ,x k )^i(x 1 )...if k (x k )dx k ...dx 1 . 

JR3 JR3 

Let us start by showing that for each k e (0, 1) and continuous function ip with compact 
support we have 

tij^lK ~ v \\c(Fx[k,x/k]) = 0, (24) 

where 

„ 5 = / ,= / *,<*,,,,)*,)*, F = {x e B? : \x\ e [k, 1/k]}. 

JR3 JR3 

Consider 

ma,£ = {£y - A)~V, m a = (£ 7 - A)"V, A ^ spec(£ 7 ). 

It is shown in [1] that converges to £ 7 in the norm resolvent sense when A ^ spec(£ 7 ). 
Thus 

\\ux, e — UxWipps) -»• as e|0, A ^ spec(£ 7 ). (25) 
Since there is a neighborhood of F (that does not include the origin) where 

(_ - A)u v = (- - A)u A = 0, 0<£<e , (26) 

when e is small enough, from standard a priori estimates for elliptic equations and (1231) 
it follows that 

|K,e - m a||c(f) — > as £ I 0, A <^ spec(£ 7 ). (27) 
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We have 



v e = f (C £ - A)-Ve A *rfA, v — — — f (£ 7 - A)-Ve A *rfA. (28) 
2™ 7r(a) 27r * 7r(a) 

Since the norm of the resolvent of an operator does not exceed the distance from the 
spectrum, we have 

II^eIIlW), ||«a||l 2 (R 3 ) < C|A| _1 , A G Q a , < 6 < E , (29) 



where fl a is the region in the complex A-plane which is to the right of the two rays starting 
at A = a and forming it/ '4 angles with the the ray [a, — oo). From here it follows that the 
contour T(a) in (|28|) can be replaced by T'(a) = dQ a : 

v £ = I u x . £ e xt d\, v = tt~t f u x e xt d\. (30) 
Zm Jr>(a) Zm Jr'(a) 

From (j2BJ) and ([29]) it follows that 

\\v\, £ \\h^f), \\u\\\h 2 (F) < c', XeT'(a), < e < e , (31) 

for some constant c', and the Sobolev embedding theorem implies that 

\\u\,e\\c(F), \\u\\\c(F) < c" , X e V (a) , 0<£<e . (32) 



This allows us to split the contour of integration in the integrals in fl30|) into two parts: 
a part where |A| >> 1 and a bounded contour. The contribution to v e and v from the 
first part can be made arbitrarily small. The Lebesgue dominated convergence theorem 
implies that the difference of the integrals over the bounded part of the contour tends to 
zero. Thus we have ( 1241) . 

It is easy to see that ff24l) implies ff23l . Let us now prove that lim £ ^ Z^(t, x ) — Z 7 (t, x). 
Let £ be an infinitely smooth function that is equal to zero for \x\ < 1 and equal to one 
for \x\ > 2. Let 



Zl — Z^' 1 + Z^' 2 , Zy — + Z^, 
where 

Z°>%x)= f %{t,x,y)Z(y)dy, Zf{t,x)=[ F y {t,x,y){l - £(y))dy, 

JR3 J-^3 

and Z , Z^ are defined similarly. Then lim £ ^ Zft' 2 (t, x) = Z (t, x) by (124"|) with ip = 
In order to estimate \Zl} x {t^ x) — Z (t,x)\, we consider the following problem: 

du 

— -C^u = h, u(p,x) = £(x), (33) 



;1 ^2 
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where h G C^°(R 3 ) depends on x only. The solution of this equation has the form 

u = Zf - -L / A _1 (£f, - X)- X he xt d\. 

^ m Jr'(a) 

Note that u(t,x) = £(x) satisfies (1331) with h = — A£/2. Thus 

The same formula is valid with Z instead of Z^ 1 and £ 7 instead of £f,. Therefore, 
\Zf(t,x)-Z\(t,x)\ = ±-\f A- 1 (^-A)- 1 (AOe A4 rfA- / \- 1 {C 1 -X)- 1 {H)e Xt d\\. 

47r ^r'(a) Jr'(a) 

Note that A£ has compact support, and therefore we can treat the right hand side in this 
formula in the same way as we treated the difference between v e and v thus obtaining 
that lim^o^' 1 ^, x) = Z 7 (t,x). 

It remains to note that the family of measures P*y, < e < 1, is tight, as follows 
from Lemma [5.11 □ 

Remark. It follows from the arguments above that for any c,t > 0, there is a con- 
stant C such that Z„(t,x) < C provided that e < c, 1/c < |x| < c and 7 < c. 

Next we prove two statements that are needed for the proofs of Theorems 12.31 and 
13.31 Note that we already demonstrated the convergence of partition functions and finite- 
dimensional distributions in Theorem 12.31 and using these facts in the proofs of Lem- 
mas ETT] and does not involve a circular argument. 

Lemma 5.1. Let t> 7 be given by |3J] ; the corresponding Hamiltonian H^ T given by (TJP 
with instead ofv, and the Gibbs measure P^'J given by (Tj|) with (3 = 1 and H^ T instead 
of Ht- For each T > and c£l, the family of measures {P^J, £ < c, 7 < c, \x\ < c} is 
tight. 

Proof. To prove tightness it is enough to demonstrate that for each r], a > there is 
< 5 < 1 such that for all u G [0, T] we have 

P*J( sup \u(s) - w(u)| > a) < 5r). (34) 

u<s<min(u+5,T) 

Let T], a > be fixed. If e > 0, the density of the measure P^J with respect to the 
measure induced by the Brownian motion starting at x is bounded from above uniformly 
in £0 < £ < c, 7 < c, 1 2 1 < c. Thus we can find < 5 < 1 such that fl34|) holds 
for Eq < £ < c. Therefore, it is sufficient to prove fl34l) under the assumption that e is 
sufficiently small so that vjL(x) = for \x\ > a/8. 
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Consider the following events in C([0, T], M 3 ), 
£j = {lu : sup |a>(s) — > a}, £j = {to : sup — u;(0)| > a/8}. 

u<s<min(u+S,T) s<min(<5,T) 

For a continuous function uj : [0, T] — > M 3 , let 

r = min(T,inf{t > : \u(t)\ > a/A}). 

Then 

P^(^) = (^(T,x))- 1 E^(exp( f v £ Mt))dt)x~ £l ) < 

Jo " 

where Ej> denotes the expectation with respect to the measure induced by the Brownian 
motion starting at the point x. Since 



E*exp(/ v e 7 (uj(t))dt)< Z £ 7 (T,x) 
Jo 



and 

K-l(Xesexp{[ T * v°(u(t))dt)) < sup E*( X£f exp( Vj (u(t))dt)), 

JO a/4<\x\<c Jo 

we have 

P;;t(^)< sup EUxejeMf v^u^dt)). 

a/A<\x\<c ' Jo 

Let 

a = min((5,inf{t > : \u(t) - w(0)| = a/8}). 
Then, since = for |x| > a/8, 

sup W T { Xe jeMf ' v*(u(t))dt)) < sup E*( X£F E^exp(/ ^H*))^)). 

a/4<|z|<c " JO a/4<|a:|<c ' JO 

Note that 

pT-a pT 

E^exp(/ ^(cu(t))rft)< sup E£exp(/ ^(o;(t))d*) = 

JO a/8<|a;|<c+a/8 JO 

sup Z*(T,x) < c(a) 

a/8<|z|<c+a/8 

for some constant c(a), where the last inequality is due to the Remark following the proof 
of Theorem 12.31 It then remains to choose 5 such that Wp{xsj) < 5i]/c(a), and estimate 
follows. □ 
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Lemma 5.2. Suppose that 7(T) is bounded and such that 7(T)vT — > +00 as T — > +00. 
Suppose that s > 0, S(T) satisfies [To]) and T is such that S(T) + s/j 2 (T) < T for 
T > T . Lei i/ie process y T {t), <t < s, be defined on the probability space C([0, T],M 3 ) 
i/ie measure P 7 (t),t by y T (t) = j(T)uj(S(T) +t/ r y 2 (T)), < t < s. Then, for each 
K > 0, the family of measures on C([0,s],M 3 ) induced by the processes y T with T > T 
and < \x\ < K is tight. 

Proof. Due to self-similarity, it is sufficient to consider the case when j(T) = 1. First 
observe that for any 77 > there is a > such that 

P X 1>T (\y T (0)\>a)< V 

for T > T and < \x\ < K. This easily follows from Lemma [3.11 in the same way that 
the convergence of finite-dimensional distributions was demonstrated in Theorem 13.31 
For a continuous function u : [0, T] — > M 3 , let 

m T (u,S) = sup \u>(ti) — u>(t 2 )\- 

|*i-*a|<5, S(T)<ti,t 2 <S(T)+s 

To prove tightness, it is sufficient to show that for each a, r\ > 0, there is 5 > such that 

P x 1T (m T (uj,S) > a) < 7] (35) 



for T >T and < |x| < K. By Theorem 12.31 

T lT {m T (u,5) >a)< liminf P^(m T (a;,5) > a). (36) 

ej.0 ' 

For a continuous function to : [0, T] — > M 3 , let 

t = mm(S(T) + s, mi{t > S(T) : \u(t)\ > a/2}). 
For < a < b < r, let 

£ r a>b = {luE C([0,r],M 3 ) : sup |w(ti) - w(t 2 )| > a}. 

|tl-*2|<5, a<tx,t 2 <b 

Then 

limmf P^(m>,5) > a) = (^(T, a;))"^ (exp^ vfaWdt)^ < 

(Z^x))-^ (exp(^ T ^(^(t))rft)E^(exp(^ T ^WMXe*-^)) < (37) 
(Z^T,x))-\Z^S(T)+s,x)) sup E*_ 5(T) (exp( / ^(u;(t))dt) X T_ sm ). 

[*[>a/2 Jo °' s 
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Since v^(x) = when x > a/4, provided that e is small enough, we have by the Markov 
property 

pT-S(T) 

sup E^„ s(T) (exp( / v £ (u(t))dt)x £ T-s(T)) < (38) 

\x\>a/2 Jo °< s 

sup 2%(T-S(T),x)P° a ({w: sup \u(h) - u(t 2 )\ > a}). 

|x|>a/4 |*i-*2|<<5i 0<t!,t 2 <s 

Recall from the proof of Theorem 12 . 31 that lim £ j Z^(t, x) = Z y (t, x). Therefore, combining 
PD, dHZD and (0HD we see that 

P~l T (m T (uj, 5) > a) < (Z J (T,x))- 1 (Z J (S(T) + s,x))x 
sup Zy(T — S(T),x)Pg({u : sup \u(ti) — u)(t 2 )\ > a}). 

|x|>a/4 \ti-t 2 \<8, 0<ti,t 2 <s 

Formula (1351) now follows from Lemma 13.11 □ 
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